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^ ■ Abstract 

'nJ" . The hyperfine structures of the ground states in the '^He'^'^ fi~ e~ and '^He^+^~e~ hehum-muonic 

atoms are investigated with the use of highly accurate variational wave functions. The differences 
O I between corresponding levels of hyperfine structure (i.e. hyperfine splittings) are determined to 

^ . very high numerical accuracy. In particular, we have found that the hyperfine structure splitting in 

ci\ 

^ i the ground state of the ^He^'^/i e atom is Au ~ 4464.55451(1) MHz, while analogous splitting 
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for the -^Re^+fi'e- atom is Au ^ 4166.43545(1). 
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I. INTRODUCTION 



In this study we report the results of highly accurate calculations of the ground 
S{L = 0)— states in the helium-muonic '^He^"*"/x~e~ and ^He^+/i~e~ atoms. Our recently 
improved methods for highly accurate variational computations allow us to construct ex- 
tremely accurate variational wave functions for these three-body helium-muonic atoms. Such 
wave functions can be used to obtain essentially exact expectation values of various bound 
state properties of these systems. In particular, the expectation values of all interparticle 
delta-functions 6ij can now be determined to very high numerical accuracy. By using these 
expectation values one can finally solve the long standing problem of hyperfine structure 
splitting in the ^He^'^//^e~ and ^He^+/i~e~ atoms 

Q-a. 

In this study we reconsider the problem of highly accurate computations of the helium- 
muonic ^B.e'^'^ fi~ e~ and ^IIe^+/i^e~ atoms. For simplicity, below we restrict ourselves to the 
analysis of the ground S{L = 0)— states in such systems. The ground state in the helium- 
muon.c atom is designated as the l»,l..-state. Analogous electton-excted states (ot l.,2.,- 
state) in the helium-muonic atoms were discussed in our earlier study [7]. Our main goal 
is to determine the hyperfine structure splitting in each of these atoms by using the results 
of most recent highly accurate computation of the expectation values of interparticle delta- 
functions. The computed values of the hyperfine structure splitting in each helium muonic 
atom must be compared with the known experimental results 8|] . For the ground states in the 
helium-muonic atoms such experiments were performed in 1970's. Analogous experimental 
observations and measurments of the hyperfine structure splitting for the electron-excited 
states in the helium-muonic atoms are difficult to perform. 

For the ground states in the ^He^"^/i~e~ and ^B.e^~^fi~e~ helium-muonic atoms the hy- 
perfine structure spllittings have been measured to very good accuracy (see references and 
discussions in jsj). The known experimental values of the hyperfine structure splitting of 
the ground states in these two atoms are 4166.41 MHz and 4464.95 MHz, respectively jsj. 
A large number of theoretical investigations of muonic-helium atoms have been conducted 
since the middle of 1970's. In many of these works the main goal was to produce highly ac- 
curate evaluation of the hyperfine structure splitting for the ground states of the '^He^+/i^e~ 
and '^He^^/i^e^ atoms. In general, to obtain the actual hyperfine structure and evaluate the 
hyperfine structure splittings one needs to detemine the expecation values of all three in- 
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terparticle delta-functions. The main computational troubles of earlier studies were related 
with the expectation values of the electron-muonic delta-function. Another complication was 
related to a slow convergence rate observed in numerical computations of the expectation 
values of the electron-nucleus delta-function. In general, it was very hard to determine the 
electron-muonic and electron-nucleus delta-functions to the accuracy better than 1 • 10^^ a.u. 
For the electron-nucleus delta-function such a 'limiting' accracy was ~ 2 ■ 10~^ a.u. Finally, 
it was almost impossible to evaluate the hyperfine structure splitting in the helium-muonic 
atoms to good numerical accuracy. A substantial progress has been achieved only in our 
work where we have applied a special Fortran pre-translator MPFUN written by David 
H.Bailey jol, 10|. By using this pre-translator we determined the hyperfine structure split- 



ting 4166.392 MHz for the ground state in the ^He^^fi^e^ atom and 4464.555 MHz for 
the ground state in the ^He^"^;Lt~e~ atom The absolute uncertainties in both cases are 
less than 10 kHz. It appears that these values were very close to the known experimental 
values of the hyperfine structure splitting measured for these atoms, which equal 4166.41 
MHz and 4464.95 MHz, respectively |8|. 

In this study we substantially improved the overall accuracy of our computations. The 
overall quality of our wave functions has also been improved. The current accuracy of the 
expectation values of interparticle delta-functions can now be evaluated as ±1 ■ 10~^ a.u. (or 
even better). Right now, we can determine the hyperfine structure splittings in both helium 
muonic atoms to much better accuracy than it was possible in 



II. THREE-BODY VARIATIONAL WAVE FUNCTION 



Our computational goal in this study is to determine the highly accurate solutions of the 
non-relativistic Schrodinger equation i/\l/(re, r^, thc) = -E^(re, 1"/^, fHe), where E < is the 
numerical parameter (energy) and H is the following Hamiltonian 
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where U, rrie and e are the reduced Planck constant, electron mass and absolute value of the 
electric charge of the electron. Everywhere below in this study, we shall use only atomic 
units, where h = l,mp = 1 and e = 1. In these units one finds 
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where masses of the muon = 206.768262 rrie and hehum nuclei MsHe = 5495.8852 rrie 
and MiHe = 7294.2996 rUe must be expressed in the electron mass mg. 

To simplify all formulas below we designate the electron, muon and helium nucleus as 
particles 1, 2 and 3, respectively. For the considered ground S{L = 0)— states in the three- 
body ^lle^~^fi~e~ and ^He^+/x^e^ muonic atoms the wave function depends upon three scalar 
interparticle coordinates, i.e. \E^(re, r^, thc) = ^('"21, ''"31, '"32)- This unknown wave function 

can be approximated to very high accuracy by using different variational expansions. 
In particular, in this work we apply the exponential variational expansion in three-body 
relative coordinates r2i,r3i and r32. For the S{L = 0)— states in the helium-muonic atoms 
this expansion takes the form 

N 

'^LM = ^Ci exp(-air2i - /^irsi - 7jr32) (3) 

i=l 

where Ci are the linear (or variational) parameters, while ai, Pi, 7^ are the non-linear param- 
eters. In general, the exponential variational expansion has significantly more complex form 
even for the S{L = 0)— states ll|. First, it is convenient to write such an expansion in the 
three perimetric coordinates ui,U2,U3 which are completely independent and each of them 
varies between and +00. In this case all non-linear parameters ai,Pi,ji must be positive. 
Second, for some three-body systems it is necessary to use the complex values for these 
parameters. It is allows one to accelerate the overall convergence rate of this variational 
expansion, e.g., for adiabatic and/or weakly-bound systems. Here we do not want to dis- 
cuss all advanteges of exponential variational expansions for three-body systems. Note only 
that the expansion, Eq.(l3]), provides a simple method HI] to construct extremely accurate 
variational wave functions for an arbitrary three-body systems (and different states in such 
systems). The overall quality of the final wave function is very good, since it can be used to 
determine all expectation values needed for analysis of the given three-body system. In par- 
ticular, such a wave function generates very accurate expectation values of all interparticle 
delta-functions which are of great interest for this study (see the next Section). 
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III. SPIN-SPIN CONTACT INTERACTION BETWEEN TWO PARTICLES 



Let us assume that we have two particles a and b which have non-zero spins and s;,, 
respectively. The non-relativistic Hamiltonian of the spin-spin interaction takes the form 

Hss = Sa-H,,= [Sa- —d^J iV 

where c is the speed of light in vacuum, n = Ms the unit vector which corresponds to the 
radius- vector r and jh is the current generated by the particle h with spin s;,. For this current 
we can write 

j, = (V X (v^*s,v^)) = (^) (n X s,) (5) 

The magnetic field generated by this spin current jf, takes the form 

Hfe = — / —Mdr i (n X (n X Si,))do = ^^\0)—Sb 6 
Sb JO ar J Sb o 

Therefore, the Hamiltonian of the spin-spin interaction is 

Hss = -^^{Sa-Sb){S{rab)) (7) 

where 5{rab) = S{ra — r?,) is the delta- function between two particles (a and b). The in- 
terparticle interaction which contains the corresponding delta-function is called the contact 
interaction. 

In atomic systems all magnetic moments are traditionally measured in Bohr magnetons 
fiB- In atomic units h = 1, rrie = 1, e = 1 the value of the Bohr magneton equals | exactly. 
In these units the formula Eq.(I7I) takes the form 

H,, = -^i^{5{rab)){Sa-Sb) (8) 

3 SaSb 

In a few-body system the total Hamiltonian of the spin-spin interaction equals to the sum 
of different two-particle Hamiltonians, Eq.dH]). In the general case, in A— body system one 
finds spin-spin terms, Eq.®, in the total Hamiltonian. In reality, many actual few- 

body systems contain identical particles and, therefore, the total number of terms in the 
Hamiltonian of spin-spin interaction is less than 
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IV. THE HYPEFINE STRUCTURE OF THE HELIUM-MUONIC ATOMS 



It follows from here that the general formula for the Hamiltonian of hyperfine structure 
{AH)h.s. of an arbitrary three-body system is written as the sum of the three following terms. 



Each of these terms is proportional to the product of the factor and expectation value 



3 

of the corresponding (interparticle) delta-funtion. The third (additional) factor contains the 
corresponding g'— factors (or hyromagnetic ratios) and scalar product of the two spin vectors. 
For instance, for the ^B.e'^'^ fi~ e~ atom this formula takes the form (in atomic units) (see, 
e.g., Q, (see also Q)) 

(AH),.,. = ^a'^^{6{rs2)){lN ■ s,) + ^«^^^(5(r3i))(I^ ■ s^) 
6 mpiTifj, 6 mprrie 

^27r^,je^ (9) 

3 mem^ 

2 

where a = |^ is the fine structure constant, and rrip are the muon and proton masses, 
respectively. The notation stands for the nucleus ^He. The factors g^,ge and are the 
corresponding (7— factors. The expression for {AH)h is, in fact, an operator in the total 
spin space which has the dimension {2sn + l)(2s^ + l)(2se + 1) = 8. In our calculations 
we have used the following numerical values for the constants and factors from Eq.(l9]): 
a = 7.2973525698 ■ 10-^ = 1836.152701me, = 206.768262me, g^ = -2.0023318414 
and ge = -2.0023193043622. The (7— factor of the helium-3 nucleus is deteremined from the 
formula: g^ = = -4.2555016, where Mn = -2.1277508 is the magnetic moments (in 
nuclear magnetons) of the helium-3 nucleus. The spin of the helium-3 nucleus is designated 
in Eq.(l9]) as In = |- 

This means that in the ground state of the helium-3 muonic atom one finds eight spin 
states. They are separated into three different groups: (1) four spin states with J = |, (2) 
two 'upper' spin states with J = |, and (3) two 'lower' spin states with J = |. Here and 
below the notation J (or J) stand for the total angular momentum J = L + S. For the 
S{L = 0)— states we have J = S, and therefore, J = S. The energy of the group of four 
states with J = | is ^ 8.29677761848-10^ MHz, while the energy of the second group of two 
'upper' spin states with J = ^ is ~ 8.29636097495-10^ MHz. The group of two 'lower' spin 
states with J = has the energy ^ —2.48899162119 - 10* MHz. The difference between 
the first and second groups of spin states is ~ 4166.43532(1) MHz. Such a difference is 
called the hyperfine structure splitting of the ground state of the helium-3 muonic atom. To 



recalculate atomic units of energy to MHz we have used the following conversion factor Ry 
= 6.579683920729-10''. 

The spin and magnetic moment of the helium-4 nucleus equals zero identically. There- 
fore, the hyperfine structure splitting in the ground state of the ^He^"''yU^e^ atom is related 
with the electron-muon spin-spin interaction. Four spin state of the electron-muon pair are 
separated in the two groups with J — 1 (three states) and J = (one state). The energy 
difference between these two groups of states is called the hyperfine structure sphtting of the 
ground state in the ^He^+//~e~ atom. With the constants mentioned above this hyperfine 
structure splitting is determined with the use of the formula 

Ai/ = 14229.178083766834 • {5{v^-^,-)) (10) 

where (5(re-^-)) is the expectation value of the electron-muon delta-function. This for- 
mula produces the hyperfine structure splittings Aiv in MHz, if the expectation value of 
delta-functions is given in atomic units. By using this formula and our expectation value 
of the electron-muon delta-function from Table II we have found that Ai/(^He^"'"/U~e~) = 
4464.55451(1) MHz. 

The uncertanties mentioned above have been determined from the corresponding uncer- 
tainties of the delta-functions. For instance, for the ^He^+z/^e" atom one finds (from Table 
II) for the (5(re-^-)) expectation value {5{r^-^,-)) ^ 3.137605340(5)-10-i a.u. Now, by us- 
ing three expectation values: 3.137605335-10-S 3.137605340-10-^ and 3.137605345-10-^ we 
have evaluated the corresponding uncertanty in the Ai/('^He^+^"e") value as < 1 - 10"^ 
MHz. 

It should be mentioned that our hyperfine structure splittings Au determined to very 
high accuracy in this study are very sensitive to variations of the basic physical constants, 
e.g., masses of particles, fine structure constant a, etc. In the case of the fine-structure 
constant a the actual changes of the numerical values of Au are quadratic upon a. This fact 
can be used to develop new experimental methods to measure the fine structure constant a 
and particle masses m^, MsHe, Mi He to very high accuracy. 
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V. CONCLUSION 



We have considered the hyperfine structures of the ground states of the hehum-3 and 
hehum-4 muonic atoms. Our predicted values of the hyperfine structure sphttings for the 
hehum-3 and hehum-4 muonic atoms are 4166.43545(1) MHz and 4464.55451(1) MHz, 
respectively. They are very close to the known experimental values 4166.41 MHz and 4464.95 
MHz, respectively jsj. Note that the corresponding experiments have been performed in 
1970's. Since then no new accurate measurments of the hyperfine splittings in these muonic 
atoms have been reported. In reality, such experiments are desperately needed to explain 
currently observed deviations between highly accurate theoretical and experimental results. 
For instance, the difference between theoretical and experimental results for the helium- 
3 muonic atom is ~ 20 times smaller than for the helium-4 atom. This has never been 
expalined. On the other hand, the new experiments are also needed to predict and explain 
some new phenomena, e.g., to determine hyperfine structure splittings in the electron excited 
states of the helium-muonic atoms. Their predicted values can be found in [3]. 
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TABLE I: The convergence of the total energies in atomic units for the ground Is^lSe— states in 
the hehum-muonic atoms. A'' is the total number of basis functions used in calculations. 

3300 -399.042 336 832 862 534 826 590 83 -402.637 263 035 135 454 018 515 17 
3500 -399.042 336 832 862 534 826 660 15 -402.637 263 035 135 454 018 585 52 
3700 -399.042 336 832 862 534 826 710 49 -402.637 263 035 135 454 018 636 94 
3840 -399.042 336 832 862 534 826 739 37 -402.637 263 035 135 454 018 666 45 
P(") -399.042 336 832 862 534 769 -402.637 263 035 135 454 004 81 

("^The best results determined in earlier calculations [8]. 

TABLE IL Convergence of the expectation values of electron-muon and electron-nucleus delta- 
functions in atomic units for the ground Is^lsg— states in the ^He^+/x~e~ and ^He^+/i~e~ helium- 
muonic atoms. The notation 1 means the electron e^, the notation 2 stands for the muon and 
notation 3 designates the nuclues of the helium-3 and helium-4, respectively. 

A(3He) (^(r2i)) (<^(r3i)) (<^(r32)) (^(r32i)) 

3300 3.1368231690-10"^ 3.206115487-10-^ 2.0149938845228-10^ 6.41400054-10^ 
3500 3.1368231739-10-^ 3.206115491-10"^ 2.0149938845221-10^ 6.41400011-10^ 
3700 3.1368231756-10-1 3.206115492-10"^ 2.0149938845218-10^ 6.41400087- 10^ 
3840 3.1368231771-10-1 3.206115494-10-^ 2.0149938845218-10^ 6.41400087- 10^ 

Nj^Re) {6{r2i)) {6{rsi)) (^(r32)) (^(r32i)) 

3300 3.1376053324-10-1 3.2063178909-10-1 2.0700137351705-10^ 6.58997578-10^ 
3500 3.1376053373-10-1 3.2063178946-10-1 2.0700137351698-10^ 6.58997637-10^ 
3700 3.1376053391-10-1 3.2063178955-10-1 2.0700137351695-10^ 6.58997612-10^ 
3840 3.1376053406-10-1 3.2063178977-10-1 2.0700137351695-10^ 6.58997614-10^ 
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